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ABSTRACT 

We present a high-frequency very long baseline interferometry (VLBl) kinematical study of the BL Lac object S5 0716-1-714 over the 
time period of September 2008 to October 2010. The aim of the study is to investigate the relation of the jet kinematics to the observed 
broadband flux variability. We find significant non-radial motions in the jet outflow of the source. In the radial direction, the highest 
measured apparent speed is ~37 c, which is exceptionally high, especially for a BL Lac object. Patterns in the jet flow reveal a roughly 
stationary feature ~0.15 mas downstream of the core. The long-term fits to the component trajectories reveal acceleration in the sub- 
mas region of the jet. The measured brightness temperature, Tg, follows a continuous trend of decline with distance, Tg oc 
which suggests a gradient in Doppler factor along the jet axis. Our analysis suggest that a moving disturbance (or a shock wave) 
from the base of the jet produces the high-energy (optical to y-ray) variations upstream of the 7 mm core, and then later causes an 
outburst in the core. Repetitive optical/y-ray flares and the curved trajectories of the associated components suggest that the shock 
front propagates along a bent trajectory or helical path. Sharper y-ray flares could be related to the passage of moving disturbances 
through the stationary feature. Our analysis suggests that the y-ray and radio emission regions have different Doppler factors. 

Key words, galaxies: active - BL Lacertae objects: individual: S5 0716-1-714 - radio continuum: galaxies - jets: galaxies - gamma- 
rays 


1. Introduction 


The combination of very long baseline interferometry (VLBI) 
images and broadband flux density variability is a unique way 
to probe the emission mechanisms near the base of jets in 
the blazar class of active galactic nuclei (AGN). The broad¬ 
band flares are often found to be connected to the ejec- 
tion of new moving emission components into the jet (e.g. , 
Marscher et all Jorstad et alll2013l : [^hinzel et al.Ll201^ 
Krichbaum et all 200lL~ and references therein). Moreover, mm- 
VLBI observations offer a unique possibility to study the struc¬ 
tural evolution in the parsec-scale region, which has been pro- 
posed to be the site of much of the high - energy emission (e.g. . 


posed to be me site ot mucn ot me nign - energy emission (e.g. , 

Runie^^, 2()j3c[ ^14t Fuhrmann et Sll2014ilMarscher et al.L 


2008H^hinzel et aLLl2012l) . Therefore, these observations have 

provided new constraints on the physical parameters of the emis¬ 
sion regions, i.e., size, brightness temperature, magnetic held, 
and motion. 


The blazar S5 0 716-H714 (z ~ 0.3, iNilsson et al.L l200l 
iDanforth et al.Ll^l2l) is a BL Lac object with a featureless opti¬ 
cal spectrum. It is one of the most intensively studied blazars be- 


cause of its extreme variability properties across the entire elec- 

troma 

snetic spectrum (e.s.. Villata e 

al.Ll2008tlFuhr 

mann et al.L 

|2()f)8l 

Rani etaLL 1201 OhIlaL 1201 3allbl 

Larionov et al.L 

2013 

). The 


broadband flux variability of the source is quite complex, with 
rapid flaring activity (on a timescale of a few hours to days) su¬ 
perimposed on top of a broad and slow variability trend on a 
timescale of ~ 1 year (iRani et al.L 1201 3at iRaiteri et'afl l2003l) . 
VLBI studies of the source show a core-dominated jet point¬ 


ing towards the north (iBach et al.L l2005l iBritzen et al.L l2009l) . 
while Very Large Array observations show a halo- l ike jet mis¬ 
aligned by ~90° on kiloparsec scales. IBritzen et alJ (l2009l) sug¬ 
gested an apparent stationarity of jet components relative to the 
core; more rec ent studies, however, have reported mo t ion as 
fast as ~40 c ( Rastorgueva et ^ 1201 It iLarionov et al.L 1201 3L 
iLister et ^ 1201 3l) . Non-radial motion and wiggling component 
trajectories have often been observed in the inner mas jet re¬ 
gion of the sourc e (IBritzen et al.Ll200^lRastorgueva et al.Ll201 It 
IRani et alll20T4l) . 


In iRani et all (l2013ah (hereafter Paper I), we presented the 
densely sampled multi-frequency observations of the source be¬ 
tween April 2007 and January 2011. These observations allowed 
us to study the broadband flaring behavior of the source and to 
probe the physical processes, location, and size of the emission 
regions. The intense optical and y-ray monitoring revealed fast 
repetitive variations (60-70 days) superimposed on a long-term 
variability trend with a time scale of ~350 days, which propa¬ 
gated down to radio wavelengths with an observed time lag of 
~65 days. A detailed investigation of the optical flares found 
a variability amplitude proportional to the flux level, which 
can be explained by a variable Doppler factor. We also found 
that the shock-in-jet model for the evolution of radio flares re¬ 
quires geometrical variations in addition to intrinsic variations 
of the source. As a possible scenario to explain the observa¬ 
tions, we suggested that the geometry significantly affects the 
long-term flux variations, which could be caused by a relativis¬ 
tic shock tracing a spiral path through the jet. To explain the 
multi-frequency behavior of an optical-y-ray outburst in 2011, 
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iLarionov et alJ (1201 3h have also suggested a shock wave propa¬ 
gating along a helical path in the blazar’s jet. 

In this paper, we use high-resolution multi-frequency VLBI 
observations to investigate the inner jet kinematics of S5 
0716-1-714, with a focus on the major radio/optical/y-ray flares 
over the time period of September 2008 to October 2010. We fo¬ 
cus on the morphological evolution of the source to investigate 
its relation to the broadband flux variations reported in Paper I, 
in particular the high-energy emission. The paper is structured 
as follows. Section 2 provides a brief description of observa¬ 
tions and data reduction. In Section 3, we report and discuss our 
results. Summary and conclusions are given in Section 4. 


2. Multi-frequency VLBI data: observations and 
data reduction 

To explore the inner jet kinematics of the source, we used the 
7 mm (43 GHz) and 3 mm (86 GHz) VLBI data obtained be¬ 
tween September 2008 and October 2010, with observations at 
26 epochs during this period. The 7 mm data are from the Boston 
University monthly monitoring program of bright y-ray blazars 
with the Very Long Baseline Array (VLBAfl The data reduc¬ 
tion was performed using standard tasks of t he Astronomical 
Imag e Processing System (AIPS) and Difmap (IShepherd et al.L 
1 19941) software. The imaging of the source (including ampli¬ 
tude and phase self-calib ration) was done using the CLEAN al- 
gorithm (H6gborr4 1974t) and SELFCAL procedures in Difmap 
dshepherd et all Il994l) . Several iterations of phase corrections 
followed by amplitude adjustments were adopted for the self- 
calibrati on process. More d etails of the data reduction can be 
found in Llorstad et al.l (l2005l) . 

The 3 mm observations were performed using the Global 
mm VLBI Array (GMVA). The GMVA is currently an array con¬ 
sisting of 14 antennas in Europe and the United States, includ¬ 
ing the 8 VLBA stations equipped with 3 mm receivers, plus the 
Effelsberg, Onsala, Metsahovi, Pico Valeta, Plateau de Bure, and 
Yebes antennas. Between October 2008 and May 2010, obser¬ 
vations were taken approximately every six months. Amplitude 
calibration is difficult at 3 mm, with atmospheric fluctuations of¬ 
ten having a major effect. We And that amplitudes are frequently 
too low after the nominal calibration in AIPS. To correct for this, 
amplitudes are scaled by comparing the flux density on the short¬ 
est baselines with single dish measurements. Details of the 3 mm 
data reduction can be found in Hodgson et al. (2014). 

We modeled the observed brightness distribution (the visi¬ 
bility amplitude and phases of the observed radio brightness of 
the source) by multiple circular Gaussian components, thus de¬ 
riving positions, flux densities, and sizes of the distinct bright 
features in the jet. The model fits were carried out for all epochs 
within the DIEMAP package, starting with a point-like model 
and fixing the position of the brightest component to (0,0). Jet 
components were added until the addition of an extra compo¬ 
nent did not lead to a significant improvement in the value of 
the fit to the uv data. The uncertainties of the model component 
parameters were determined by comparing the parameter ranges 
obtained after performing model fits with a different number of 
model components. We used at least four different model fits for 
each epoch to obtain the parameter uncertainties. However, the 
uncertainties depend also on the self-calibration and data editing, 
and on the brightness and size of the components. In addition, 
the uncertainties increase with increasing distance from the core. 
This is accounted for by following an independent approach for 


VLBA-BU-BLAZARS, http://www.bu.edu/blazars 


the error estimation proposed bv iKrichbaum et ^ (Il998h . The 
fitted model parameters for all of the epochs are listed in Table 
lAl where we only list the formal errors obtained from the first 
method. Eigures[T] and |2] show the fitted circular Gaussian com¬ 
ponents superposed on the clean images. 

To investigate the kinematics in the jet of S5 0716-H714, 
the individual model components were identified following the 
assumption that the changes of the flux density, distance from 
the VLBI core, position angle, and size should be small for the 
time period between adjacent epochs. In order to prevent a po¬ 
tentially large systematic error arising from the incorrect cross- 
identification of moving features from epoch to epoch, the sim¬ 
plest scheme was adopted while identifying the jet-features. A 
self-consistent cross-identification is proposed using all avail¬ 
able model-fit parameters. These cross-identifications are not 
necessarily unique, especially if the source evolves in a more 
complex manner than assumed here. 


3. Results and discussion 

In this section we present the flux and spatial evolution of the 
bright radio emission in the jet of S5 0716-1-714. Special atten¬ 
tion is given to the kinematics of the inner jet (< 2 mas from the 
core). We also investigate possible correlations between the jet 
kinematics and the broadband flux variations reported in Paper 
1 . 

3.1. Jet kinematics 

Eigures [1] and |2] display a sequence of VLBI maps convolved 
with a natural beam. The model fit parameters are given in Table 
lAl In the following sub-sections, we will discuss in detail the ap¬ 
parent motions of components and brightness temperature gra¬ 
dient in the parsec-scale jet. 

3.1.1. Component motion 

Eor the kinematic study of individual components, we choose 
the VLBI core as a reference point, fixed to coordinates (0,0). 
However, it is important to note that the absolute position of 
the reference point i.e. core may change e.g. due to changes 
in opacity or instabilities in the jet; but, the results presented 
here are not sensitive to that as we are interested in the relative 
motion of components w.r.t. the core. The VLBI data were 
fitted by using different numbers of components during different 
activity stages over the time span considered here. During this 
period, we identified a total of 12 components - Cl, Kl, Al, 
A2, A3, and NCI to NC7 - in addition to the core, CO (see 
Fig- El)- Eigure[3]plots the evolution of the distance of different 
knots from the core, and Eig. |4] shows their trajectories in the 
XY-plane projected on the sky. 

Stationary Features: 

The sequence of VLBI images allow us to investigate the 
component motion along the jet as a function of time. Most 
of the components exhibit significant motion down the jet, 
except Al, A2, and A3. The latter are comparatively stable in 
position, although they exhibit significant scatter in position 
angle exceeding the uncertainties (see Eig. |4]i. Nevertheless, 
their radial distance from the core remains at ~0.15 mas. 
Statio nary features are a common chara c teristi c in AGN jets 
(e.g. lEromm et all l2013al: iJorstad etakl 120011 : iBritzen et all 
l2010h . In straight jets, these can be produced by recollimation 
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Fig. 1. 43 GHz VLBA images (contours) and Gaussian components of S5 0716-1-714 at different epochs between September 2008 
and October 2010. The date of observation is given at the bottom of each map (shifted Julian day, JD-2454000). The fitted model 
parameters for each epoch are given in Table The cross at the bottom of each map represents the FWHM of the restoring beam 
(0.25x0.18 mas). Contour levels correspond to -0.003 0.003 0.006 0.012 0.024 0.048 0.096 0.19 0.38 0.76 1.54 3.07 6.14 12.28 
24.58 49.15 98.30 mJy/beam. 


shocks, instabilities [magnetohydrodynamic (MH D), and/or 
Kelvin-Helmholtz (KH)1, or magnetic p inches (see iMarschetl 
I2009I: iHardeel l20ok iMeier et all l200l[ for details). Bends in 
the jet can also cause quasi-stationary features, either because 
the jet turns more into the line of sight, thus increasing the 


Doppler beaming facto r, or due to t h e form ation of a shock 
that deflects the flow (lAlberdi et al.L Il993l) . Numerical sim¬ 
ulations indicate that when a moving knot passes through a 
standing re-collimation shock, the components blend into a 
single feature, then split up after the collision with no lasting 


3 
























B. Rani et al.: Jet kinematics and high-energy emission connection in 0716-1-714 



Relative R. A. (mas) 


Fig. 2. The inner jet structure of S5 0716-1-714 from GMVA observations at 3 mm. The cross at the bottom of each map represents 
the FWHM of the restoring beam (0.07x0.04 mas). Contour levels correspond to -0.004 0.004 0.008 0.016 0.032 0.064 0.12 0.25 
0.51 1.02 2.04 4.09 8.19 16.38 32.77 65.54 mJy/beam. 



Fig. 3. Evolution of core separation of the individual components 
as a function of time. 


changes in the prop e r moti on of the moving knot (iGomez et al.L 
Il997t iFromm et akl l2012l) . During the interaction, the feature 
may move a short distance downstream before returning to its 
previous position. Therefore, it is quite possible that Al, A2, 
and A3 are the same standing feature, which becomes disturbed 
by moving features. We find that the kinematics of three new 
components NC2, NC3, and NC6, are significantly different 
before and after the interaction. The components move more 
slowly (fiapp < 10c) as they approach the stationary feature, 
then faster (fiapp ^ 20c) after the interaction. Such behavior is 
con sistent with the stati onary feature representing a bend in the 
jet (lAlberdi et al.L 1 19^ . Indeed, a bend at ~0.15 mas is evident 
in the 3 mm maps of the source (see Fig. |2l). This suggests that 
the observed stationary features are standing oblique shocks 
that cause the flow of the jet to bend at ~0.15 mas owing to a 
transverse pressure gradient or an impact with an interstellar 
cloud. 

Moving Features: 

The other components. Cl, Kl, and NCI to NC7, exhibited sig- 



X (mas) 

Fig. 4. Trajectories of the jet knots in the X-Y plane, where X = 
r.cos(0 -t- 90) and Y = r.sin(0 + 90) (r is the radial separation from 
the core and 0 is the position angle with respect to an imaginary 
north-south line drawn through the map center). 


nificant motion in the radial direction. Component Cl represents 
very faint emission at a core separation of >1 mas. It seems to 
be the remnants of a component ejected in earlier epochs. The 
bright feature Kl c ould be a trailing co mponent that forms be¬ 
hind a strong shock (lAgudo et al.L 1200 H) as it is identified at four 
epochs after the ejection of a new component (NCI). 

As seen in Fig.lH all components follow curved trajectories 
in the XY-plane projected on the sky. In addition, the trajecto¬ 
ries differ from one component to the next. The wiggling trajec¬ 
tories in the jet might be a signature of hel i cal rn otion. Helical 
jet models (e.g.. lGomez et al.pl994lHardeeLl2006l) represent the 
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jet as an inhomogeneous flow of plasma interacting with ambient 
matter. Magnetohydrodynamic or Kel vin-Helmholtz instabilities 
(iHardeel l2006t iPerucho et al.L l2006h can explain the bending 
and helical structures, and possibly also less regular wiggles ob¬ 
served in parsec-scale radio jets. 

In order to investigate component motions in the jet of S5 
0716H-714, we use the XY positions of the components relative 
to the core to fit the trajectories with polynomials of different 
order. The resulting fitted radial separation of each component, 
r(f) = yjxitf + is shown in Fig. |5] We find that a linear 
function is sufficient to fit the trajectories of components NCI, 
NC4, NC5, and NC7. However, a second-order polynomial is 
required to fit the trajectories of components NC2, NC3, and 
NC6, which corresponds to an apparent acceleration. Two mul¬ 
tiple straight lines needed to fit different portions of the trajecto¬ 
ries of these components in Fig.|5]provides a clear demonstration 
of the inadequacy of simple linear fits. The dashed curves repre¬ 
sent the quadratic fits. 

For each component, the fits yield an average proper motion 
and a mean speed. Using the estimated angular speed fi, we have 
computed the kinematic parameters of the jet, e.g., the apparent 
speed, Papp, and Doppler factor, 5. We derive the apparent speed 
of the components, ySapp, from the angular speed, /r, using 


A 


app 


^^dL 

C(1 H-z)’ 


( 1 ) 


where dp is the luminosity distance and z the redshift of the 
source. The luminosity distance dp corresponding to z=0.3 is 
dp — 1600 Mpc for a ACDM cos mology with Q,„ - 0 .27, Qa - 
0.73, and Hq-1\ km s ' Mpc ' (ISnergel et al.L 120031) . The an¬ 
gular and apparent speed of the individual components are given 
in Table [T] The calculated Papp values range from 6c to 37c, the 
latter of which is unusually high for a blazar, especially for a 
BL Lac object. 

For the non-ballistic components, NC2, NC3, and NC6, the 
quadratic fits represent significant acceleration along the jet 
axis. Formally, we obtain p - 4.90+1.55 mas yr“^ for NC2, p 
- 2.45+0.80 mas yr“^ for NC3, and p = 0.94+0.36 mas yr“^ 
for NC6. Such apparent acceleration is a commo nly observed 
featur e of many objects in the MOJAVE sample (iLister et al.L 
1201 3h . for which a statistically significant tendency for acceler¬ 
ation in the base of jets has been found. A standing shock slows 
the flow while deflecting it, after which the flow can accelerate 
back to the velocity it had upstream of the shock (iGomez et al.L 
Il997h . However, since the apparent speed was relatively slow 
in the upstream region, this seems not to fit the case of S5 
0716+714. The data are more consistent with the scenario 
that, after a moving component passes through the stationary 
feature, it bends towards or away from our line-of-sight. It then 
appears to move faster because its velocity vector subtends 
an angle closer to that which maximizes apparent motion, 
arccos(;S), where p is the velocity of the component divided 
by the speed of light. An alternative possibility is that the 
underlying flow of the jet systematically accelerates outward. 
Theoretical models involve strong magnetic fields associated 
with the putative supermassive black hole/accretion disk system 
that play a k ey role in the initia l accele r ation and collimation 
of Ae jet ( Blan dford & PavneL Il982t iBlandford & ZnaieS. 
Il977t iMeier et ~ . l200lh . The conversion of Poynting flux to 
flow energy is gradual and may persist out to parsec scales 
(e.g.. lSikora et al.Ll2005h . Under this explanation, however, it is 
difficult to explain the straight-line kinematic evolution of some 
of the components. 



Fig. 5. Evolution of radial separation of different components 
from the core. The solid lines represent fits with linear functions, 
while the dashed curves are quadratic fits. 


Variability Doppler factor: 

A varia bility Doppler factor using VLBI jet kinematics was de¬ 
fined bv iJorstad et al.l (l2005h as 

dvLBI - —T- (2) 

where s is the angular size of the component, defined as 1.6x0 
for a Gaussian with EWHM diameter 0 measured at the epoch 
of maximum brightness, and dp is the luminosity distance. The 
timescale of variabilit y is defined as At^ar — dtjlniS maxIS min) 
(iBurbidge et al.LIl974l) . where dt is the time separation in years 
between maximum (Smax) and minimum (Smm) flux densities. 
The estimated values of 6vpbi for the individual components, 
listed in Table [1] range from 6 to 21. Using the broadband 
flux density and spectral variability study of the source over the 
same time period (Paper I, Sect. 3.4), we obtain a range of self- 
consistent values for the Doppler factor. The different indepen¬ 
dent approaches suggest that 6 >20, which is consistent with 
upper range of the values of dvpBi obtained here. 


3.1.2. Ejection Epochs 

Back-extrapolation of the components’ motion allows us to es¬ 
timate the time of zero separation from the core (i.e., the “ejec¬ 
tion” time. To). As shown in the previous section, the compo¬ 
nents follow curved trajectories; therefore, we back-extrapolate 
the component trajectories in the XY-plane to determine the 
epoch of (0,0) separation. The error estimates are not very 
straightforward to calculate for components following non- 
ballistic trajectories. We use the estimated error on the fitted pa¬ 
rameters for a given function to calculate its uncertainty. Back- 
extrapolation of the two envelopes of the function provide us 
with the error in Tq. The calculated ejection times for the indi¬ 
vidual components are listed in Table [T] 


3.1.3. Brightness temperature gradient in the jet 

The redshift-corrected brightness temperature (Tg_„ij) of the 
bright emiss ion features can be approximated using the follow¬ 
ing relation (iJorstad et al.L 120051) : 


TB,obs = 1.22 X 10 


12 


' comp 


(1 + z) dlmnp V- 


K, 


(3) 


where S comp is the component flux density in Jy, dcomp is the 
EWHM of the circular Gaussian component in mas, and v is the 
observing frequency in GHz; for the calculations here, we have 
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Table 1. Physical parameters of the radio emission 


Component 

yu (mas/yr) 

Pappic) 

^var 

To (days) 

NCI 

1.30+0.04 

24.88+0.88 

7.03+0.35 

783™ 

NC2, 

0.45+0.00 

8.71+0.13 

13.84+0.69 

900![“ 

NC2„ 

1.93+0.14 

37.06+2.75 



NC3, 

0.61+0.03 

11.69+0.64 

6.02+0.30 

933!]® 

NC3„ 

1.87+0.17 

35.78+3.29 



NC4 

0.87+0.04 

16.81+0.95 

6.91+0.34 

1047! 

NC5 

0.85+0.03 

16.30+0.59 

11.72+0.58 

1079!^ 

NC6, 

0.34+0.04 

6.53+0.77 

15.33+0.76 

1092!2s 

NC6„ 

0.99+0.04 

19.05+0.92 



NC7 

0.48+0.04 

9.20+0.84 

20.96+1.04 

1203!]' 


I and u denote, respectively, lower and upper values of the component 

speeds. 



Fig. 6. Component brightness temperature as a function of ra¬ 
dial separation from the core. The solid line represents the fitted 
power law with index e = 2.43 + 0.36. 


only used 43 GHz data. The calculated brightness temperatures 
{Ts.obs) of the components generally decline with radial sepa¬ 
ration from the core (see Fig. This decaying trend can be 
approximated by a power law, Ts^obs ~ fj/, with e = 2.43+0.36, 
where r/e, is the radial distance along the jet axis. The solid line 
in Fig. [^represents the fitted power law. 

We place the observed brightnes s temperature gradient i n 
the context of the shock-in-jet model dMarscher & GeapIlOSSh . 
In the common picture of this model, a relativistic shock 
propagates down a conical jet, slowly expanding adiabatically 
while maintaining shock conditions during propagation. While 
the shock propagates down the jet, it undergoes three ma¬ 
jor evolutionary stages dominated by inverse Compton, syn¬ 
chrotron, and adiabatic energy losses. As a result, the observed 
brightness temperature declines as a power law, Ts.obs ^7^ 

jei 

The value of e can be derived from the spectral evolution 
of the radio emission (iLobanov & Zenstii. 119991 : iFromm et al.L 
I2013bl). Following Tb tx and S(v) oc 

(iLobanov & Zensu^ll999l) . we obtain 

n-s-a(s+l) b(s + 3) ^ 

^Compton — Q ^ 1 ,^,) 


^synchrotron ~ 




1)] ^ b(s + 3) ^ 


(5) 


and 

2(5-2.)-3fl(.+ l) X^ + 3) ,, 

^adiabatic ~ ^ 

O 1 

where a, b, and s parametrize the variations along the jet axis of 
the Doppler factor, 6 oc magnetic field, B oc and power 
law distribution of energy of the emitting electrons, N(y) oc 
Assuming a constant Doppler factor, i.e., b -0, and using a typ¬ 
ical value of i = 2 (corresponding to optically thin synchrotron 
spectral index = -0.5), and a = 1 to 2 (for toroidal and poloidal 
magnetic fields, respectively), we obtain ecompton - 1.3 to 1.6, 
^synchrotron ~ 4.2 tO 5.7, and €odiabatic — 5.2 tO 4.7. 

The observed brightness temperature gradient in S5 
0716+714 (Fig. |8]l has a slope equal to 2.43+0.36. The slope 
of the observed intensity gradient rules out the simple assump¬ 
tions of a constant Doppler factor. We therefore consider a vari¬ 
able Doppler factor along the jet axis. We find that the observed 
value of e is consistent with the adiabatic loss phase if 5 oc 
with Z? ~ 0.12 to 0.40, which corresponds to a moderate variation 
in the Doppler factor. However, with relatively larger variations 
in the Doppler factor b ~ 0.47 to 0.74, the observed value of e 
would agree with synchrotron loss phase as well. 

The evolution of the flux density and frequency (5„, and v,„) 
at the spectral turnover from synchrotron self-absorption of the 
radio flares, discussed in Paper I (Section 3.3.2), also suggests a 
variation of 6 along the jet axis. In that study, we found that 6 
changes as r* with b - 0.7 during the rise and b - 2.6 during 
the decay of the first radio flare. The evolution of the second flare 
was governed hy b = 0.4 during the rising phase and b - -2.0 
during the decay (see Paper I for details). Therefore, it is evi¬ 
dent that the two flares require a different dependence of 6 on 
distance from the core, which is also inconsistent with what we 
obtained for the brightness temperature gradient (b = 0.12 to 
0.74). This can be explained in two ways. First, the single-dish 
radio flux correlates well with the VLBI core flux (see Section 
l3.2l for details). This implies that the radio flares originate within 
the core. The core is more compact and farther upstream than the 
jet components, with a comparatively higher brightness temper¬ 
ature, and can follow a different dependence of 6 on rjet- The 
individual components trace the intensity gradient further down¬ 
stream of the jet, and therefore can individually have a different 
6 - rjet dependence. A second possibility is that the brightness 
temperature gradient plot corresponds to an averaged behavior 
of the two radio flares. 

3.2. Jet kinematics and broadband emission 

For the comparison of the observed broadband flaring behav¬ 
ior with the jet kinematics, we have used the flux density light 
curves from Paper I, where the details of observations and data 
reduction are given. Broadband flux light curves are shown in 
Fig. [7] Panel (a) presents the weekly averaged y-ray photon flux 
light curve. Panel (b) shows the X-ray light curve, while panel 
(c) displays the optical V-band light curve. Panel (d) presents the 
radio (single-dish) and 43 GHz cor^l flux density curves plotted 
in the same panel. Panel (e) exhibit flux density light curves of 
the individual components shown in different colors, with the 
same color used to shade their respective ejection times. The 
major outburst between JD’ ~840 and 1150 is marked with a 
horizontal arrow in panel (b). The repetitive optical flares with a 

^ We used a mean value of the m easured optically thin spectral in¬ 
dex, = 0.4 iRani et all 1201 3ali to scale the 86 GHz flux density 
measurements to those at 43 GHz. 
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Fig. 7. Broadband flux light curves of S5 0716-1-714: (a) weekly averaged y-ray flux light curve at E>100MeV, (b) X-ray light curve 
at 0.3-10 keV, (c) optical V-band light curve, (d) single-dish (43 and 230 GHz) and 7 mm VLBI core flux light curves, and (e) flux 
density curves of moving jet components shown in different colors; the same color is used to shade their respective ejection times, 
i.e., epochs of zero separation from the core. The green dashed lines indicate the average times when moving components cross the 
stationary feature at ~0.15 mas. The major outburst is marked with a horizontal arrow; the faster repetitive optical flares are labeled 
as “0” to “9” and sharp y-rays flares as “A” to “J” (see text for details). 


typical duration of ~120-140 days are label as “0” to “9”, and 
the broadband flares are denoted “G” for y-rays, “X” for X-rays, 
“O” for optical and “R” for radio followed by the number ad¬ 
jacent to them; e.g., we have flares GO to G9 at y-rays, X6 at 
X-rays, GO to 09 at optical, and R6 and R9 at radio frequen¬ 


cies. The sharp y-ray flares are labeled as “A” to “J”, and have a 
duration of ~10-20 days. 

Although VLBI observations of S5 0716H-714 are missing at 
the peak of the two major radio flares (JD’ ~ 1150-1200 and 
~1300), there is still an apparent one-to-one correlation between 
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Table 2. Component ejection epochs and broadband flares 


Component 

Flares observed at the time of 


zero separation 

component crosses 


from the core 

stationary feature 

NCI 

OO/GO 

GOB 

NC2 

02/G2 

G4C 

NC3 

03/G3 

04D 

NC4 

05E/G5E 

06G/G6G 

NC5 

05E/G5F 

G6I 

NC6 

06/G6 

- 

NC7 

R6 

09J/G9J 


the core flux and single-dish radio flux measurements. This sug¬ 
gests that the single-dish radio flares are produced in the core, 
which agrees with the high mean core-to-jet flux ratio, which 
ranges from 1.3 to 21.4 with a mean of 4.9+4.0. The broadband 
flux cross-correlations (Paper I) also indicate that the flux varia¬ 
tions at y-ray (and optical) and single-dish radio frequencies are 
correlated such that the peak of the major optical/y-ray outburst 
leads the peak of the radio outburst by a two-month time period. 

3.2.1. Component ejections and their connection with 
broadband flares 

Because of the superposition of different modes of flaring ac¬ 
tivity, the broadband flaring behavior is quite complex in the 
source. The major outburst (marked with a horizontal arrow in 
Fig.0 is accompanied by faster repetitive optical/y-ray flares 
(labeled as “0” to “9”) and sharp y-rays flares (labeled as “A” 
to “J”). The jet kinematic study during this period suggests sig¬ 
nificant changes in the jet morphology of the source; in total, 
seven new components were ejected. It is therefore plausible 
that the ejection of new components could be connected to the 
broadband flares. We have identified a stationary feature at a dis¬ 
tance of ~0.15 mas from the core (Section 3.1.1). Interaction of 
a moving component with the stationary feature could amplify 
the magnetic h eld and accelerate particles to produ ce a multi¬ 
waveband flare (lMarscheiil2014l : lFromm et al.Ll2012h . 

In Fig. [T] the shaded areas represent the ejection times of 
new jet components and the green dashed lines mark the average 
epochs when a moving component crosses the stationary feature 
at ~0.15 mas. In Table |2l we list the new components (col. 1), 
along with the flares that coincide with their respective ejection 
times (col. 2). Flares observed closer in time to the passage of 
a moving component across the stationary feature are listed in 
col. 3. We note that most of the sharp y-ray flares (B, C, D, G, 
I, and J) are observed when a moving component crosses the 
stationary feature, while most of the repetitive faster flares (0, 2, 
3, 5, and 6) coincide with the ejections of new components from 
the core. 

In Paper I, we found that the high-energy (optical and y-ray) 
emission is significantly correlated with the radio flares such that 
the major optical/y-ray outburst propagates down to radio fre¬ 
quencies (R6 flare) with a time delay of ~65 days, following a 
power-law dependence on frequency with a slope ~0.3. We also 
noticed that the evolution o f radio flare R6 is in a ccordance with 
a generalized shock model (IValtaoia et al.L[l99^ . This suggests 
that the origin of the major outburst could be related to a distur¬ 
bance/shock that moves down the jet. 

We interpret the broadband flaring events in the follow¬ 
ing manner (see Fig. A disturbance at the base of the jet 
could cause a shock wave to form and propagate down the jet 


(IValtaoia et al.L 119921 : iMarscher & Geaii Il985h . which could be 
related to the onset of the major optical/y-ray outburst (ID’ 
~840) and disturbs the jet outflow, resulting in the ejection of 
two new components (NC2 and NC3, Fig. [8^). The core is op¬ 
tically thick at that point, as the shock forms upstream of the 
core; for this reason, we do not see much activity in the core. 
Repetitive optical/y-ray flares and the curved trajectories of the 
associated components favor the propagation of the shock along 
a bent trajectory, possibly a helical path. The moving shock 
reaches its maximum strength at ID’ ~1150 (Fig. [8]b), close 
to the peak of the major optical/y-ray outburst (see Section 3.3, 
Paper I for details), which also coincides with the onset of radio 
flare R6. At this point, the shock reaches the core, which starts 
to become optically thin. After this, further disturbances pass¬ 
ing through the core correspond to the ejections of three new 
components (NC3, NC4, and NC5). Major radio outburst R6 is 
observed ~65 days later, accompanied by the ejection of com¬ 
ponent NC7 (Fig.[8]c). As the moving features cross the the sta¬ 
tionary feature at ~0.15 mas, at first, the stationary feature gets 
destructed after the interaction, and later it appears again. The 
times of passage of moving components through the stationary 
feature appear to coincide with the sharp y-ray flares in most 
cases. Therefore, the sharp y-ray flares seem to be produced via 
the interaction of moving features with the stationary one. 

X-ray flare X6 coincides with a similar increase in the flux 
of the 43 GHz core. Optical flare 07 peaks about a month later, 
which roughly coincides with the expected passage through 
the stationary feature of the same disturbance that caused flare 
X6. However, no bright moving knot is apparent in the next 
VLB A image ~ 20 days after 07. Also, the dramatic decrease 
in the millimeter-wavelength flux during this interval implies 
that the disturbance is strongly quenched following flare 07. 
There are two anomalies in this sequence of events: there is 
no optical flare as the disturbance crosses the core, and there 
is no y-ray activity during the entire time span covering flares 
X6 and 07. The first indicates that, for reasons that are un¬ 
clear ( e.g., non-favorable magnetic held orientation relative to a 
shock: ISummerlin & Baringll2012h . electrons with energies high 
enough to produce optical and y-ray emission were not produced 
in the core region at this time. The second implies that there was 
a dearth of seed photons for inverse-Compton scattering ~0.15 
mas from the core during flare 07. 


3.2.2. Doppler beaming of radio and gamma-ray emission 
regions 

The observed time lag between y-ray and 43 GHz core flares 
can be used to estimate the separation b etween the two emissio n 
regions by using the following equation (iFuhrmann et alil2014l) : 




/^app^ 

sin0 


(7) 


where 0 is the viewing angle of the source, jSapp is the appar¬ 
ent jet speed, and Ty r is the o bserved time lag. Using the ob¬ 
served value Tyy = 82+32 days (IRani et al.Ll2014li . we have Aryy 
= (O.O69±O.O26)y6app(sin0)“' pc. 

The measured apparent speeds in the source cover a wide 
range. The current observations suggest slower motion (<8 c) on 
sub-mas scales (~0.3 mas), which later get accelerated to speeds 
as fast as 40 c on mas scales. The source also has a prominent 
bent jet structure at sub-mas scales (see Fig.|2]i; as a result, the 
sub-mas and mas regions could have different viewing angles. 
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Black 

hole 


Stationary 

feature 



High-energy 
Emission region 


7 mm VLBI core 


(a): On set of the major outburst (JD' = 840) 


Black 

hole 


7 mm VLBI core 


(b): The major optical/gamma-ray outburst (JD' = 1150) 



Stationary 

feature 


Stationary 

feature 



7 mm VLBI core 


(c) : The major radio outburst (JD' = 1190) 

Fig. 8. Sketch for the proposed scenario in the BL Lac S5 
0716+714 (not to scale). Note that north is to the right in this 
sketch. A shock propagates down the jet along a helical path, 
producing a major optical/y-ray outburst upstream of the 7 mm 
VLBI core and later a major radio outburst accompanied by ejec¬ 
tion of components from the core. 


For ySapp = 6-8 c, the viewing angl^E G, is 6-9°; this gives ^ 
= 2.9-4.4 pc. Using ySapp = 40 c, we have 6 - 1.4° and Aryy 
~ 110 pc. E ven for a larger viewin g angle, 0=10°, we obtain 
Ary r ~16 PC. iPushkarev et aP (1201 2l) estimated that the 15 GHz 
core is located at a distance of 6.68 pc from the jet base/apex. 
Using the technique to align the optically-thin emission based on 
model-fitted components and morphological similarities, we de¬ 
rive a shift of 0.042+0.013 mas (0.19+0.06 parsec) between the 
15 GHz and 43 GHz core, which puts the location of the 43 GHz 
core at a distance of ~6.5 pc from the jet-apex. Therefore, the 
distance between y-ray emission region and 43 GHz core cannot 
be larger than 6.5 pc. For a typical viewing angle 9 » 5°, Ary^ < 
6.5 pc gives ySapp < 8 c. However, we note that for a larger view¬ 
ing angle {9 ~ 10°), /lapp could be as high as 10-12c but not ~ 40c 
unless the viewing angle is larger than 30°, which is unlikely to 
be the case for S5 0716+714. The results therefore suggest that 
the high-energy and radio emission regions have different appar¬ 
ent speeds and probably different viewing angles as well. This 
implies that the y-ray and radio emission regions have different 
Doppler factors. 

0 = 5i«^‘(l/j8app) 


4. Summary and conclusion 


The monthly sampled mm-VLBI monitoring of S5 0716+714 
allows us to study the inner jet kinematics and to investigate its 
relation to the observed broadband flux variability. Our study re¬ 
veals significant non-radial motions in the jet outflow, including 
variations in the orientation of the sub-parsec-scale jet (see Rani 
et al. 2014 for details) and possibly helical trajectories of moving 
components. 

In radial directions, the individual components exhibit ex¬ 
treme apparent speeds as high as ~37 c. An apparent speed 
of ~20 c seems to be more t ypical in S5 0716+714 on parsec 
scales. Recently. iLister et"^ (120131) reported an apparent speed 
of 43.6+1.3 c. These values would even be extreme for a quasar. 
Therefore, moving components in the jets of BL Lac objects 
can be just as fast as in quasars, despite flux-limited surveys at 
lower resolution (e.g. iLister et SI l2009l) finding that the mean 
apparent velocity of BL Lac objects is lower than for quasars. 
Further mm-VLBI studies are needed to determine the extent 
to which angular resolution and the variety of spectral classes of 
BL Lac objects (low- to high-frequency peaks of the synchrotron 
spectral energy distributions) skew the apparent velocity distri¬ 
bution. It is possible that the fast BL Lacs are actually quasars in 
disguise, with their broad emission lines swamped by the non- 
thermal optical continuum. On the other hand, emission lines 
have never been observed in S5 0716+714. 

The bright components in 0716+714 follow a continuous 
decay trend in the observed brightness temperature (Tb o^^ oc 
which implies a variation in Doppler factor along 
the jet axis. The spectral evolution of radio flares (see Paper I for 
details) also suggests variations in Doppler factor. A significant 
correlation between y-r ay flux and varia tions in orientation of 
the sub-parsec scale jet (iRani et al.Ll2014l) supports the possibil¬ 
ity that the observed flux variations in S5 0716+714 are strongly 
influenced by Doppler factor variations, and therefore may have 
a geometrical origin. 

Our study suggests a connection between jet kinematics and 
the observed broadband flaring activity. We find that the origin 
of major optical/y-ray/radio outbursts in the source seems to be 
related to a disturbance, such as shock wave, propagating down 
the jet. The repetitive, faster optical/y-ray flares, superposed on 
the major outbursts, and non-ballistic motions of the associated 
jet components can be interpreted via helical jet models (as dis¬ 
cussed by, e.g., Blandford & Pavnel 1982t Blandford & Znaiefl 


I 1977 I : ICamenzind & Krockgnto^ . 1992 ) and/or magnetohy 
drodynamic instabilities ( Hardeel 20061) . Our analysis favors that 
the major optical/y-ray flares in S5 0716+714 are produced up¬ 
stream of the 7 mm VLBI core. The sharp y-ray flares, how¬ 
ever, seem to be produced when the moving components cross a 
stationary feature 0.15 mas from the core. Finally, our analysis 
suggests that the upstream y-ray emission region has a different 
Doppler factor than observed in the parsec-scale radio jet. 
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Table A.l. Results from Gaussian Model fitting and component parameters 


Epoch (JD') 
Date 

^ peak 

(Jy/beam) 

r 

(mas) 

e 

n 

0 

(mas) 

Comp" 

Wavelength 

720 

1.670±0.084 

0±0 

0±0 

0.044±0.002 

core 

7 mm 

10.09.2008 

0.139±0.007 

0.098±0.005 

39.90±2.00 

0.034±0.002 

Al 



0.023±0.001 

0.714±0.036 

13.28±0.66 

0.427±0.021 

Cl 


787 

1.670±0.084 

0±0 

0±0 

0.038±0.002 

core 

7 mm 

16.11.2008 

0.188±0.009 

0.103±0.005 

29.29±1.46 

0.071 ±0.004 

Al 



0.022±0.001 

1.516±0.076 

12.33±0.62 

0.487±0.024 

Cl 


822 

2.044±0.102 

0±0 

0±0 

0.034±0.002 

core 7 mm 


21.12.2008 

0.082±0.004 

0.139±0.007 

26.83±1.34 

0.101±0.005 

NCI 



0.013±0.001 

1.700±0.085 

13.57±0.68 

0.538±0.027 

Cl 


856 

2.167±0.108 

0±0 

0±0 

0.034±0.002 

core 7 mm 


24.01.2009 

0.405±0.020 

0.108±0.005 

8.06±0.40 

0.045±0.002 

Al 



0.052±0.003 

0.292±0.015 

15.92±0.80 

0.171±0.009 

NCI 



0.025±0.001 

1.572±0.079 

11.01±0.55 

0.879±0.044 

Cl 


885 

1.717±0.086 

0±0 

0±0 

0.035±0.002 

core 

7 mm 

22.02.2009 

0.420±0.021 

0.125±0.006 

6.64±0.33 

0.051 ±0.003 

Al 



0.053±0.003 

0.383±0.019 

9.89±0.49 

0.201±0.010 

NCI 



0.016±0.001 

1.873±0.094 

13.47±0.67 

0.642±0.032 

Cl 


924 

1.038±0.052 

0±0 

0±0 

0.018±0.001 

core 

7 mm 

02.04.2009 

0.552±0.028 

0.036±0.002 

56.28±2.81 

0.036±0.002 

NC2 



0.187±0.009 

0.150±0.008 

13.26±0.66 

0.062±0.003 

Al 



0.051 ±0.003 

0.460±0.023 

7.60±0.38 

0.163±0.008 

NCI 



0.020±0.001 

1.672±0.084 

8.84±0.44 

0.692±0.035 

Cl 


963 

0.354±0.018 

0±0 

0±0 

0.011 ±0.001 

core 

3 mm 

11.05.2009 

0.104±0.005 

0.028±0.001 

73.33±3.67 

0.012±0.001 

NC3 



0.043±0.002 

0.084±0.004 

47.52±2.38 

0.025±0.001 

NC2 



0.048±0.002 

0.159±0.008 

54.57±2.73 

0.030±0.001 

Al 



0.030±0.001 

0.441 ±0.022 

16.93±0.85 

0.047±0.002 

K1 



0.019±0.001 

0.628±0.031 

39.92±2.00 

0.044±0.002 

NCI 


982 

0.482±0.024 

0±0 

0±0 

0.038±0.002 

core 

7 mm 

30.05.2009 

0.118±0.006 

0.109±0.005 

29.31±1.47 

0.074±0.004 

NC2 



0.035±0.002 

0.525±0.026 

5.61±0.28 

0.192±0.010 

K1 


1004 

0.548±0.027 

0±0 

0±0 

0.039±0.002 

core 

7 mm 

21.06.2009 

0.075±0.004 

0.104±0.005 

31.10±1.56 

0.097±0.005 

NC3 



0.022±0.001 

0.594±0.030 

8.26±0.41 

0.148±0.007 

K1 



0.564±0.028 

0±0 

0±0 

0.035±0.002 

core 

7 mm 

1040 

0.166±0.008 

0.080±0.004 

39.50±1.98 

0.040±0.002 

A2 


27.07.2009 

0.038±0.002 

0.165±0.008 

19.10±0.96 

0.067±0.003 

NC3 



0.010±0.001 

0.690±0.035 

11.10±0.56 

0.187±0.009 

K1 


1060 

0.663±0.033 

0±0 

0±0 

0.039±0.002 

core 

7 mm 

16.08.2009 

0.125±0.006 

0.079±0.004 

56.40±2.82 

0.040±0.002 

A2 



0.030±0.002 

0.188±0.009 

23.50±1.18 

0.066±0.003 

NC3 



0.019±0.001 

0.494±0.025 

6.77±0.34 

0.181±0.009 

NC2 



0.006±0.000 

1.098±0.055 

2.52±0.13 

0.104±0.005 

NCI 


1091 

0.873±0.044 

0±0 

0±0 

0.045±0.002 

core 

7 mm 

16.09.2009 

0.172±0.009 

0.112±0.006 

48.70±2.44 

0.044±0.002 

NC4 



0.041 ±0.002 

0.388±0.019 

11.50±0.58 

0.199±0.010 

NC3 



0.015±0.001 

0.764±0.038 

9.97±0.50 

0.147±0.007 

NC2 


1117 

0.356±0.018 

0±0 

0±0 

0.012±0.001 

core 

3 mm 

12.10.2009 

0.110±0.005 

0.029±0.001 

73.71±3.69 

0.015±0.001 

NC6 



0.041 ±0.002 

0.084±0.004 

48.90±2.45 

0.017±0.001 

NC5 



0.042±0.002 

0.158±0.008 

53.49±2.67 

0.024±0.001 

NC4 



0.021±0.001 

0.502±0.025 

17.16±0.86 

0.040±0.002 

NC3 


1119 

0.835±0.042 

0±0 

0±0 

0.039±0.002 

core 

7 mm 

14.10.2009 

0.122±0.006 

0.100±0.005 

47.00±2.35 

0.041±0.002 

NC5 



0.025±0.001 

0.563±0.028 

11.50±0.58 

0.197±0.010 

NC3 



0.013±0.001 

0.875±0.044 

15.00±0.75 

0.189±0.009 

NC2 


1121 

1.142±0.057 

0±0 

0±0 

0.040±0.002 

core 

7 mm 

16.10.2009 

0.175±0.009 

0.100±0.005 

48.10±2.41 

0.041±0.002 

NC5 



0.038±0.002 

0.607±0.030 

10.30±0.52 

0.248±0.012 

NC3 



0.007±0.000 

1.201 ±0.060 

20.30±1.02 

0.352±0.018 

NCI 


1125 

1.204±0.060 

0±0 

0±0 

0.039±0.002 

core 

7 mm 

20.10.2009 

0.243±0.012 

0.103±0.005 

43.20±2.16 

0.052±0.003 

NC5 



0.036±0.002 

0.648±0.032 

9.79±0.49 

0.238±0.012 

NC3 



0.010±0.001 

1.177±0.059 

9.19±0.46 

0.391±0.020 

NCI 



JD' = JD - 2454000, 

Speak '■ The integrated flux in the component, 

r : The radial distance of the component center from the center of the map, 
6 : The position angle of the center of the component, and 
(j ): The FWHM of the component. 
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Table 1 continued. 


Epoch (JD') 
Date 

^peak 

(Jy/beam) 

r 

(mas) 

e 

n 

<P 

(mas) 

Comp" 

Wavelength 

1131 

1.446±0.072 

0±0 

0±0 

0.035±0.002 

core 

7 mm 

25.10.2009 

0.262±0.013 

0.115±0.006 

40.70±2.04 

0.052±0.003 

NC5 



0.029±0.001 

0.675±0.034 

10.90±0.55 

0.224±0.011 

NC3 



0.013±0.001 

1.256±0.063 

14.30±0.72 

0.227±0.011 

NCI 


1152 

4.144±0.207 

0±0 

0±0 

0.032±0.002 

core 

7 mm 

16.11.2009 

0.397±0.020 

0.075±0.004 

31.60±1.58 

0.064±0.003 

NC6 



0.077±0.004 

0.254±0.013 

27.30±1.37 

0.099±0.005 

NC4 



0.018±0.001 

0.736±0.037 

10.70±0.54 

0.214±0.011 

NC3 


1207 

1.880±0.094 

0±0 

0±0 

0.043±0.002 

core 

7 mm 

10.01.2010 

0.206±0.010 

0.108±0.005 

20.30±1.02 

0.060±0.003 

NC6 



0.025±0.001 

0.380±0.019 

17.20±0.86 

0.224±0.011 

NC4 



0.015±0.001 

0.817±0.041 

9.89±0.49 

0.271±0.014 

NC3 


1239 

1.503±0.075 

0±0 

0±0 

0.035±0.002 

core 

7 mm 

11.02.2010 

0.219±0.011 

0.046±0.002 

48.90±2.45 

0.041±0.002 

NC7 



0.064±0.003 

0.138±0.007 

15.50±0.78 

0.072±0.004 

NC6 



0.031 ±0.002 

0.355±0.018 

11.90±0.60 

0.280±0.014 

NC5 



0.007±0.000 

0.951 ±0.048 

12.07±0.60 

0.152±0.008 

NC3 


1262 

1.453±0.073 

0±0 

0±0 

0.035±0.002 

core 

7 mm 

06.032010 

0.477±0.024 

0.082±0.004 

38.90±1.94 

0.041±0.002 

NC7 



0.050±0.002 

0.198±0.010 

6.09±0.30 

0.072±0.004 

NC6 



0.016±0.001 

0.550±0.027 

12.60±0.63 

0.254±0.013 

NC4 



0.011±0.001 

1.464±0.073 

7.88±0.39 

0.374±0.019 

NC3 


1327 

0.902±0.045 

0±0 

0±0 

0.020±0.001 

core 

3 mm 

10.05.2010 

0.194±0.010 

0.065±0.003 

71.62±3.58 

0.032±0.002 

A3 



0.019±0.001 

0.560±0.028 

14.54±0.73 

0.024±0.001 

NC7 



0.213±0.011 

0.134±0.007 

70.29±3.51 

0.063±0.003 

X 



0.023±0.001 

0.200±0.010 

45.00±2.25 

0.021±0.001 

NC5 


1336 

1.543±0.077 

0±0 

0±0 

0.041±0.002 

core 

7 mm 

19.05.2010 

0.359±0.018 

0.109±0.005 

58.30±2.92 

0.062±0.003 

A3 



0.043±0.002 

0.435±0.022 

12.40±0.62 

0.232±0.012 

NC6 



0.009±0.000 

1.624±0.081 

15.20±0.76 

0.323±0.016 

X 


1362 

0.783±0.039 

0.000±0.000 

0±0 

0.041±0.002 

core 

7 mm 

14.06.2010 

0.261±0.013 

0.110±0.005 

49.50±2.48 

0.075±0.004 

A3 



0.020±0.001 

0.687±0.034 

12.00±0.60 

0.251±0.013 

NC6 



0.005±0.000 

1.781±0.089 

15.60±0.78 

0.171 ±0.009 

X 


1410 

1.389±0.069 

0.000±0.000 

0±0 

0.040±0.002 

core 

7 mm 

01.08.2010 

0.438±0.022 

0.070±0.003 

70.20±3.51 

0.053±0.003 

A3 



0.080±0.004 

0.250±0.013 

29.70±1.49 

0.173±0.009 

NC7 



0.014±0.001 

0.620±0.031 

8.66±0.43 

0.181±0.009 

NC6 


1458 

2.974±0.149 

0±0 

0±0 

0.033±0.002 

core 

7 mm 

18.09.2010 

0.535±0.027 

0.081 ±0.004 

73.60±3.68 

0.040±0.002 

A3 



0.030±0.002 

0.185±0.009 

50.50±2.53 

0.050±0.003 

X 



0.022±0.001 

0.324±0.016 

29.80±1.49 

0.176±0.009 

NC7 



0.032±0.002 

0.693±0.035 

10.90±0.55 

0.235±0.012 

NC6 



a : Identification of the individual components. If a component appeared only in a single epoch, we labeled it with X. 
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